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■ Abstract. We prove the reality of the perturbed eigenvalues of some PF symmetric 

Hamiltonians of physical interest by means of stability methods. In particular we study 
' 2-dimensional generalized harmonic oscillators with polynomial perturbation and the one- 

dimensional x^{ixY for — 1 < e < 0. 
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1 Introduction and statement of the results 



A basic mathematical question in the framework of VT symmetric quantum mechanics concerns 
\^ ' the reality of the spectrum of the considered Hamiltonian. One crucial issue is to formulate con- 

• ditions for the reality of (part of) the spectrum of Hamiltonians depending on a real parameter e 

i of the form 

Hie) = Ho + eV, e G R. (1) 

o ■ 

In [H [2] results have been obtained in the case when the unperturbed Hamiltonian Hq is self- 
adjoint with discrete spectrum, and in [3] in the more general case when Hq is VT symmetric but 
^ \ not necessarily selfadjoint, with discrete spectrum. The common framework of these papers is 

^ ■ perturbation theory, which allows to obtain the reality result only for small values of the coupling 



constant e. However this is not a limitation: in many cases as e varies, critical phenomena occur, 
such as a spontaneous breaking of the VF symmetry and the appearance of complex eigenvalues 
caused by the crossing of energy levels of Hq. The most elementary example of this phenomenon 
is represented by a 2 x 2 matrix of the form 

H.ie)=('' (2) 

with 61,62 G M, recently examined in [H [5]. -f^i(e) is VT symmetric, i.e. Hi{e){VT) = 
(VT)IIi{e), if the parity operator V is the unitary involution defined by 



V 



1 
-1 
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and T is the antilinear involution of complex conjugation. The eigenvalues of Hi{e) are 
A(e) = l(ei + 62) ± iV(ei-e2)2-4e2 

and they are real if and only if |e| < ^^^^^^ . The matrix ([2]) can be rewritten in the perturbative 
form ([T]) 

so that it can be put in the framework of Theorem 1.2 of |T], whose applicability goes far beyond 
the case of finite dimensional operators (it includes in fact classes of Schrodinger operators with 
bounded perturbation potentials, also in dimension greater than one): this theorem ensures the 
reality of the spectrum of Hi{e) for |e| < where denotes the norm of the bounded 

perturbation operator W, which in the present case is 1 since W = ( ^ n I • "^^^ conclusion 



is that the perturbation argument of [1] yields actually the reality of the spectrum for all the 
allowed values of e. 

In a similar way one can analyze VT symmetric matrices of the form 

with e, 6 S M, which are VT symmetric with 
\l 

For a physical interpretation of these matrices see [HIS] for applications to optics, and Section U 
of the present paper for an application to classical mechanics. Complex Hamiltonians of type ([5]) 
are a particular case of those considered in [6]. In turn the Hamiltonians ^ are a particular 
case of those examined in [7]. 

The results obtained in [Tl[2] deal with the case when the perturbation potential V is bounded 
and, under suitable assumptions on Hq and V, they guarantee the reality of the entire spectrum 
of H{e) for je| < (5/|[H^|[, where 5 = | inf |A„ — Ami and (t{Hq) = {A„ : n G N} is the (discrete) 

spectrum of Hq. In some cases, not only the reality of the eigenvalues can be proved, but also 
the similarity to a selfadjoint operator [SI [9]. In [3] the authors analyze the case when V is 
unbounded (and in general not even relatively bounded with respect to Hq), obtaining a weaker 
result which guarantees that the perturbed eigenvalues of H(e), converging to those of Hq as 
e — 7> 0, are real for |e| sufficiently small. In |10] this result has been extended to classes of 
Hamiltonians of the form 

H{e) = H + iW„ e > 

acting in L^(M), where H = denotes the operator associated with the one dimensional 

harmonic oscillator and € C''(M) is an odd real- valued function: Wf_{—x) = —We{x) € M, 
Vx S M, Ve > 0. Such generalization of the result of [3] is based on the fact that in the 
dependence on the perturbation parameter e is the most general one and not only of linear type 
as in ([1]). The aim of this paper is to show how simple extensions of the results of [3l [10] allow 
us to treat models of considerable interest in the context of VT symmetric quantum mechanics. 
More precisely, extending the results of [3] to the case of Schrodinger operators in dimension 
greater than 1 it is possible to examine models of the following type 



H2ie) = --j^ - + oofxl + uj^xl + iex[x2, x = {xi,X2) eM? , (4) 
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where we assume r, s G N, r + s odd, and OJ 1,002 > 0. This type of Hamiltonians, which 
represent a natural generalization of the quantum Henon-Heiles model, has been studied also 
in [T2] . We will analyze the Hamiltonian Q in Section [2] and we will prove the reality of 
the eigenvalues Ani,n2(£) of H2{e) generated by the unperturbed ones, i.e. by the eigenvalues 
Ani,n2(0) = (2ni + l)u}i + (2n2 + l)uj2, Vni,n2 S No := {0,1,2,...}, such that A^^^nje) ^> 
Ani,n2(0) as e ^ 0, in the case of non-resonant frequencies wi, uj2- 

Then, generalizing also the results of |10] we will examine the basic model 



The Hamiltonian ([5]) has been examined and its physical interest discussed by several authors 
[131 [m [151 [16] and the reality of its spectrum for e > has been proved in [l7j (see also [H]); for 
e < numerical results indicate the appearance of complex eigenvalues (for a recent discussion 
see e.g. [H]), in the sense that, as |e| increases while e moves from to —1, the perturbed 
eigenvalues of H^{e), generated by the unperturbed ones A„(0) = 2n + 1, Vn E No, become 
complex. In Section [3] we will prove that as long as |e| stays sufficiently small, those eigenvalues 
remain real. We will actually treat the case of small |e| independently of the sign of e: for e > 
we recover the known result proved in |17j . 

Remark 1. As in [31 [10], the results of Sections [2] and [3] are obtained by proving the stability 
of the unperturbed eigenvalues. For the notion of stability see e.g. [191 [20]. Iii particular, let us 
recall that if {H{e) : |e| < eo} is a family of closed VT symmetric operators in a Hilbert space % 
and if a simple (i.e., non degenerate) eigenvalue A of -ff(O) is stable with respect to the family 
{H{e) : |e| < eo}, then for |e| sufficiently small there exists one and only one eigenvalue A(e) 
of H{e) near A such that 

A(e) A, as e 0. 

Now, recalling that the eigenvalues of a VT symmetric operator come in pairs of complex 
conjugate values, we can assert that the uniqueness of A(e) implies its reality: in fact, if A(e) is 
not real then there are two distinct eigenvalues, A(e) and A(e), and not just one as stated above. 

Assume now that %= L'^{W^), d>\, and that H{e), \e\ < eo, is formally given by 



H{e) =p^ + Vi + W,, 

d 

where = —A = - ^ x = (xi, . . . , Xk) G M'^ and Vi, G L^^^^R'^). Moreover we assume 



that H{e) has discrete spectrum, i.e. each point in the spectrum is an isolated eigenvalue with 
finite multiplicity (for a review of the notion of multiplicity of an eigenvalue see also [3]). Let 
H*{e) denote the adjoint of i^(e) and assume that both H{e) and H*{e) have C^{R'^) as a core. 

As recalled in [10| the stability of the eigenvalues of H{0) in the framework of the stabil- 
ity theory developed by Hunziker and Vock in [21] holds if the following three conditions are 
satisfied. 




e 



-1< e < 1. 



(5) 



fc=i 



(1) For allueC§°{R'^) 



lim 



H{e)u-H{0)u\\ =0 



(6) 



and 



lim 



H*{e)u- H*{0)u\\ =0. 



(7) 
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(2) There exist constants a,b > and 7, I7I < tt/2, such that 

{u,p'^u) < a{cosj?R.{u, H{e)u) + sin j'is{u, H{e)u) + b{u,u)} (8) 

for all u e C^{R'^). 

(3) For any z € C, there exist positive constants 5, no and ei such that 

dn{z,Hie))>5>0, (9) 
for all n > uq and \e\ < e\, where dn{z,H{e)) := dist(z, A/'n(e)) and 

Afnie) := {{u, H{e)u) : u G D{H{e)), \\u\\ = 1, u{x) = for |x| < n} 

is the so called "numerical range at infinity" {see [3]). 

Here D{H{e)) denotes the domain of H{e) and {u,v) := J^au{x)v{x)dx denotes the standard 
scalar product in L'^{W^). For future reference we summarize this result in the following criterion, 
proved in [21] (see also [TO]). 

Proposition 1. The eigenvalues of H{0) are stable w.r.t. the family H{e), \e\ < eo, if the above 
conditions (l)-(3) are satisfied. 

In Sections [2] and [3] we will prove that the operator families H2{e) and H^{e) satisfy this 
criterion and, in view of Remark [1] above, this is enough to ensure the reality of the perturbed 
eigenvalues of H2{€) and H^{e). 

2 The operator family H2{e): stability and reality 
of the perturbed eigenvalues 

The Schrodinger operator H2{e), e G R, introduced in Section[T]is defined as the closure in L^(M^) 
of the minimal operator defined on Cq°(M'^), we denote its domain D(H2{e)), and its adjoint 
is given by H2{—e). It has discrete spectrum (see e.g. [20], Theorem XIII. 64; we note that the 
extension to the present non-selfadjoint case is straightforward) and the unperturbed eigenvalues, 
i.e. the eigenvalues of -^2(0), are given by 

-^ni,n2(0) = (2ni + l)a;i + (2n2 + l)a;2, Vni,n2 G Nq. 

We assume that the frequencies oji,oj2 > are non-resonant, i.e. the equation kiuji + A:2a;2 = 0, 
ki,k2 G Z, is satisfied if and only if ki,k2 = 0. Then each eigenvalue Ani,n2(0) is simple, i.e. 
the corresponding eigenspace has dimension 1. Moreover we assume that r, s G N are not both 
even. Then H2{e) is VT symmetric, where T is complex conjugation and the parity operator V 
is defined as follows: if r and s are both odd V changes the sign of one coordinate, i.e. it 
can be either {Viu){xi,X2) = u{—xi,X2) or {V2u){xi,X2) = u{xi, —X2), Vx = {xi,X2) G M^, 
Vti G L^(M^). If r and s are not both odd a suitable choice for V is {V3u){x) = u{—x), Vx G M^, 
Vii G L^M.'^). 

We now prove that conditions (l)-(3) are satisfied. 

Proposition 2. Conditions (l)-(3) are satisfied by the operator family H2{e), e G M. 

Proof. (1) Vn G Cq°(M^) let K = supp(ti) denote the support of u; since K is compact in R^, 
there exists a constant M > such that < M, Vx = (xi, X2) G K. Therefore 

\\H2{e)u - H2{^)uf = \\H:^{e)u - Hl{{))u\\^ = [ |kx^x|u(x)p(ix 

Jk 

<e'^M^ [ \uix)\'^dx = e'^M^\\uf, Ve G R, 
Jk 

and this yields 1^ and ([7]). 
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(2) Vn G C^(M2) we have 

^{u,H2{e)u) = {u^p^u) + ijjI{u,x\u) +a;|(n, x^n) > {u,p'^u). 

Thus, ([8|) is satisfied with 7 = 0, a = l,6 = 0. 

(3) Let z G C and uj := mm{uji,uj2}- Let u € D{H2{e)) be such that ||n|| = 1 and u{x) = 
for |a;| < n, n N. Then 

\z- {u,H2ie)u)\ > \{u,H2ie)u)\ - \z\ > ^{u,H2{e)u) - \z\ 

\ 2/ 2 \ , 2/ 2 \ I I \ 2/ I |2 \ I I \ 2 2 i i 
> UJi [U, XiU) + UJ2{U, X2U) — \z\>UJ [U, \x\ U) — \z\ > LO n — \z\. 

Thus, dist(z, A/'n(e)) > u'^n- \z\, Ve € M, Vn € N, whence lim dn{z,H{e)) = +00, Ve € M and 
this proves ([9]). ■ 

Corollary 1. Near each unperturbed eigenvalue A„^^„2(0); "'i)'^2 G I^Oi of H2{0) there exists 
one and only one eigenvalue A„^^„2(^) of H2{e) for e G M, |e| small, converging to A„^^„2(0) as 
e ^ 0. Moreover A„^^„2(e) is real. 

Proof. The first statement follows from Proposition [T] (applicable after Proposition [2]), the 
definition of stability of eigenvalues and the fact that each eigenvalue Ani,n2(0) of H2{0) is 
simple. The reality of Ani,n2(e) follows from its uniqueness, as anticipated in Remark [H ■ 

3 The operator family H3{e): stability and reality 
of the perturbed eigenvalues 

We first consider the minimal operator in L^(M) defined on C^(M) by the formal expression ([5]) 
for |el < 1, i.e. Vn G C^(M) we set 



H3{e)u = -u" + x^{ixfu. (10) 
Then (jlOp can be rewritten as follows 

H3{e)u = -u" + e^^'S'^^)*^^ \x\^+'u = -u" + cos (f e) \x\^+'u + (sign x)i sin (f e) jx^+^n, 

f 1, if X > 
where sign x = < ^ . „ 

[ -1, if X < 0. 

Then H3{e) is closable and the domain of its closure (still denoted H^^e)) is D{H3{e)) = 
i?^(IR) n L'(|xp^'^), Ve: |e| < 1. Moreover the closed operator i^3(e) has compact resolvent and 
therefore discrete spectrum, and it is VT symmetric if T is, once again, the complex conjugation 
operator and V is the parity operator defined by {Vu){x) = n(— x), Vu G L"^' 



The unperturbed operator H3{0) = + x^, where p'^ = — corresponds to the one- 
dimensional harmonic oscillator and its eigenvalues A„(0) = 2n + 1, n G Nq, are simple. Let 
us now proceed in analogy with Section [2] and prove that the operator family H^i^e) satisfies 
conditions (l)-(3). 

Proposition 3. Conditions (l)-(3) are satisfied by the operator family H-i{e), —1 < e < 1. 
Proof. (1) Vn G C^(M.) let K = supp(n). Then 

\\H3{e)u - H3{0)u\\^ = [ {e^''^'''=^''^\x\'^+' - x'^Y\u{x)\^dx. (11) 

Jk 

Since K is compact there exists a constant c > such that the integrand expression in the 
right hand side of (fTTj) can be bounded from above by c|u(x)p, Vx G K, Ve G]— 1, 1[. Then, 
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by Lebesgue's dominated convergence theorem the r.h.s. of (|lip converges to zero as e — ?■ and 
this proves ([U]). A similar argument proves d?]). 

(2) Vn G C^(R2) we have 

^{u, H^{e)u) = {u^p^u) +cos(§e) {u,\x\'^^^u) > {u,p^u). 

Thus, as for the case of H2{e), ([8]) is satisfied with 7, 6 = and a = 1. 

(3) Again, in analogy with the argument used for H2{e), let 2; € C and u G D{H^{e)) be such 
that = 1 and u{x) = for \x\ < n, n £ 'N. Then 

\z - {u, Hs{e)u) \ > ?ft{u, H3{e)u) - \z\ > cos (f e) {u,\x\'^+'u) - \z\ > ^ - \z\ 

for |e| < ei := |. ■ 

Now, with an argument analogous to that used to prove Corollary [1] we obtain the following 
result. 

Corollary 2. Near each eigenvalue A„(0) = 2n + 1, n € No, of H^iQi) there exists one and only 
one eigenvalue Xn{^) of H^^e) for e €] — |e| small, converging to An(0) as 0. Moreover 
A.„(e) is real. 

4 Conclusions 

Concerning Section [1] we remark that 2x2 matrices can characterize complexified classical 
mechanical systems (see e.g. [22]) related to Schrodinger operators dH) for r = s = 1, i.e. 

P1+P2+ U^fxi + U!2^2 + i£XlX2. 

Indeed the classical equation of motion can be written as 

where Xk, k = 1,2, denotes the second time derivative of Xk- It is straightforward to study the 
eigenvalue problem for the 2x2 matrix in (jl2p . which corresponds to ([2]) with = 2w|, k = 1, 2. 
The eigenvalues are 

A±(e) = {ul + ul) ± ^{u:l-u^lf-e\ 

We see that the condition |e| < — i^H, which yields the reality of the spectrum in quantum 
mechanics, also yields the reality of the frequencies of the normal oscillation modes in classical 
mechanics. The quantum eigenvalues are given by 

A„i,n2(e) = (2ni + l)A+(e) + (2^2 + l)A_(e), Vni,n2 € Nq. 

Expanding the functions A-i-(e) in powers of e we obtain the Rayleigh-Schrodinger pertur- 
bation expansion (RSPE) for the eigenvalues A„j^„2(e). All expansions clearly converge for 
|e| < — wll, therefore the radius of convergence of the RSPE coincides with the threshold of 
transition between real and complex spectrum. 

In a similar way one proves that the complexified classical Hamiltonian corresponding to the 
2x2 matrix ([3]) for e, 6 € M 

Pi+P2 + \{^ + ie)x\ + \{e- ie)x1 + hxiX2 
admits real normal modes with real frequencies if |e| < \b\. 
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As far as Section [2] is concerned we remark that a suitable generalization regards polynomial 
perturbations of the d-dimensional harmonic oscillator, d > 2. In addition, still in dimension 
d = 2 one could try to study the case of resonant frequencies, not with stability methods (which 
fail because of the degeneracy of the unperturbed eigenvalues), but proving the reality of the 
Rayleigh-Schrodinger perturbation expansion and its summability. 

Finally, concerning Section[3]an open question is the rigorous proof of the existence of complex 
eigenvalues for —1 < e < 0, supporting the existing numerical results, and the analysis of the 
case e < —1. 
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